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WELL-POSEDNESS FOR A MULTI-DIMENSIONAL VISCOUS 
LIQUID-GAS TWO-PHASE FLOW MODEL 

CHENGCHUN HAO AND HAI-LIANG LI 

Abstract. The Cauchy problem of a multi-dimensional {d ^ 2) compressible vis- 
cous liquid-gas two-phase flow model is concerned in this paper. We investigate the 
global existence and uniqueness of the strong solution for the initial data close to 
a stable equilibrium and the local in time existence and uniqueness of the solution 
with general initial data in the framework of Besov spaces. A continuation criterion 
is also obtained for the local solution. 



1. Introduction 



The models of two-phase or multi-phase flows have a very broad applications of 
p I ■ hydrodynamics in industry, engineering, biomedicine and so on, where the fluids un- 

-^ . der investigation contain more than one component. Indeed, it has been estimated 

(-H ! that over half of anything which is produced in a modern industrial society depends, 

to some extent, on a multi-phase flow process for their optimum design and safe op- 
erations. In nature, there is a variety of different multi-phase flow phenomena, such 
as sediment transport, geysers, volcanic eruptions, clouds and rain. In addition, the 
-^ , models of multi-phase flows also naturally appear in many contexts within biology, 

>- 1 ranging from tumor biology and anticancer therapies, development biology and plant 

^ ■ physiology, etc. The principles of single-phase flow fluid dynamics and heat transfer 

are relatively well understood, however, the thermo-fluid dynamics of two-phase flows 
is an order of magnitude more complicated subject than that of the single-phase flow 
^ I due to the existence of moving and deformable interface and its interactions with two 

phases [3, 14, 15]. 

We consider the drift-flux model of two-phase flows in the present paper, which is 
principally developed by Zuber and Findlay (1965), Walhs (1969) and Ishii (1977). 
The basic idea about drift-flux models is that both phases are well mixed, but the 
relative motion between the phases is governed by a particular subset of the flow 
5^ ! parameters. In the case of liquid-gas fluids, it relates the liquid-gas velocity difference 

to the drift-flux (or "drift velocity" ) of the vapor relative to the liquid due to buoyancy 
effects. In general, the drift-flux models consist of two mass conservation equations 
corresponding to each of the two phases, and one equation for the conservation of 
the mixture momentum, and are particularly useful in the analysis of sedimentation, 
fluidization (batch, cocurrent and countercurrent), and so on ([13, 18, 23, 24]). 

The Cauchy problem to a simplified version of the viscous compressible liquid-gas 
two-phase flow model of drift-flux type in M'^ [d ^ 2), where the gas phase has not 
been taken into account in the momentum equation except the pressure term and the 
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equal velocity of the liquid and gas flows has been assumed, reads 
frit + div(mu) = 0, 

fit + div(nu) = 0, (1.1) 

(rfiu)t + div(mu ® u) + 'S/P{fh, n) = flAu + (/i + A) Vdivu, 

with the initial data 

(m,n,u)|i=o = ("^o,^o,Uo)(x), in M"^, (1.2) 

where fh = aipi and n = Q.gPg denote the liquid mass and the gas mass, respectively. 
The unknowns a;, a^ G [0, 1] denote the liquid and gas volume fractions, satisfying the 
fundamental relation: ai + Ug = 1. The unknown variables pi and pg denote the liquid 
and gas densities, satisfying the equations of states pi = pifi + {P — Pifl)/af, pg = P/a^, 
where ai and ag denote the sonic speeds of the liquid and the gas, respectively, and 
Pifi and pifi are the reference pressure and density given as constants, u denotes the 
mixed velocity of the liquid and the gas, and P is the common pressure for both 
phases, which satisfies 

P{m,n) = Co (—b{m,h) + \/h'^{m, h) + c{fri, h) j , (1.3) 

with Co = af/2, k^ = pi^o — Pifi/af > 0, oq = ct^/^f and 

6(m, n) = ko — fh — aofi, c{fh, h) = Akoaofi. 

p and A are the viscosity constants, satisfying 

p>0, 2p + d~X^0. (1.4) 

For the one-dimensional case, the existence and/or uniqueness of the global weak 
solution to the free boundary value problem was studied in [9, 11, 21, 22] where the 
liquid is incompressible and the gas is polytropic, and in [10] where both of two fluids 
are compressible. However, there are few results for mult i- dimensional cases except 
for some computational results [17]. As a generalization of the results in [10], the 
existence of the global solution to the 2D model was obtained in [20] for small initial 
energies. In [19], a blow-up criterion for the 2D model was proved in terms of the 
upper bound of the liquid mass for the strong solution in a smooth bounded domain. 

One of the main results of the present paper is the existence and uniqueness of 
the global strong solution to the Cauchy problem (1.1)-(1.2) under the framework of 
Besov spaces, for all multi-dimensions d ^ 2, provided that the initial data are close 
to a constant equilibrium state. The other result is the local well-posedness and the 
continuation criterion to the Cauchy problem with general initial data. Because of the 
similarity of the viscous liquid-gas two-phase flow model to the compressible Navier- 
Stokes equations, we can apply some ideas adopted in the proof of well-posedness for 
the compressible Navier-Stokes equations to deal with the two-phase flow model. It is 
Danchin who first makes important progress in applying the Littlewood-Paley theory 
and Besov spaces to sovle the existence and uniqueness for the compressible Navier- 
Stokes equations or barotropic viscous fluid in [5, 8] and for the flows of compressible 
viscous and heat-conductive gases in [6, 7]. However, it is non-trivial to apply di- 
rectly the ideas used in single-phase models into the two-phase models because the 
momentum equation is given only for the mixture and that the pressure involves the 
masses of two phases in a nonlinear way, which makes it rather difficult to obtain the 
estimates of the masses and the mixed velocity (m, n, u) in Lebesgue spaces L^ with 
respect to the time. In addition, it seems impossible to get the estimates of fh and n 
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from the system simultaneously due to the strong coupling among the corresponding 
terms, even if we change the variables (m, n) linearly. 

To overcome these difficulties in global well-posedness theory, we make use of a 
nonlinear variable transform so as to separate the two mass variables from each other, 
which enable us to decompose the original system into a transport equation and 
a coupled hyperbolic-parabolic system. To be more precisely, we first divide the 
momentum equation by m (which supplies additional information) and take a new 
variable n = ao{h/m — n/m) for some constants n and rh. This makes the resulted 
equation for n a homogeneous transport equation with the velocity u, and the expected 
estimates of the new variable depend only on the mixed velocity. Then, we remove 
the linear terms involving n from the momentum equation so as to separate linearly 
the equation about n from the others, which can be done by virtue of the variables 
changes with a careful choice of coefficient factors. Finally, to establish the a priori 
estimates for the global existence theory, we deal with the linearized system directly 
instead of separating the velocity into the compressible and incompressible parts. 

As for the local well-posedness theory for general data, we need to reformulate the 
original system and deal with the resulted nonlinear system directly, and in terms 
of the improvement of the a priori estimates on the densities, we can generalize the 
local well-posedness result in [1, 8] to the two-phase flow model (1.1) with the specifled 
pressure function. 

Before stating the main results, we introduce some notations. Throughout the 
paper, C (or c) stands for a harmless constant, and we sometimes use A < B to 
stand for A ^ CB. B^ and 5*'* denote usual homogeneous Besov spaces and hybrid 
Besov spaces, respectively; L°°(i?*'*) and C{B'^'^) are mixed time-spatial spaces, see 
the appendix for details. Let us now introduce the functional spaces which appear in 
the theorems. 

Definition 1.1. For T > and s G M, we denote 
E^= {(m,n,u) ■.neC{[0,T];B'-^''{R'^)) 

m e C([0, T]; B'-^'^iR"^)) n L\[0, T]; B'+^'^iR"^)) 
uG {C{[0,T]; B'-^R'')) n L\[0,T]; B'+^R''))^}, 

and 

||(m,n, u)||£;j = ip||x,ooQo,T];B=-l.=) + ll"^ll L°°([0,T];B=-i.=) + l|u|lz,°°{[0,T];_B^-i) 

+ ||''^||li([0,T];B=+i.^) + ||u||li([o,T];B''+1)- 

We use the notation E'^ ii T = +oo, changing [0, T] into [0, oo) in the deflnition 
above. 

Definition 1.2. Let a G [0, 1] and T > 0, denote 
F^:=(C([0,T];5'^/2,d/2+a))i+i 

X (C([0, T]; 5'i/2-i,rf/2-i+a) n ^iQo, T]; 5'i/2+l,d/2+l+a))rf_ 

Now, we state the global well-posedness results briefly as follows. For more infor- 
mation about the solution, one can see Theorem 2.1 in the second section. 

Theorem 1.1 (Global well-posedness for small data). Let d ^ 2, n '^ 0, rh > 

(1 — sgnn)/co; /i > and 2fl + dX ^ 0, in addition, fl + X > if d = 2. There exist two 
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positive constants a and Q such that if rho — rh, Uq — n E i^'^/2-i.'i/2 ^^^c? uq G B'^^'^^^ 
satisfying 

||mo — rfi\\Qd/2-i.d/2 + llfio — n\\Qd/2-i.d/2 + ||uo||^d/2-i ^ a, (1-5) 

then the following results hold 

(i) Existence: The system (1.1) has a solution (m, n, u) satisfying 

rfi-m,fi-neC (M+; B''I^~^^'"^) , u G C (M+; B'"^~^) , 
and moreover, 

II (a(m — m) + bao{h/m — n/fh),h/ih — n/fh, u)||£;d/2 

^(5(||mo — m||5d/2-i,d/2 + ||fio — n||^d/2-i,d/2 + ||uo||5d/2-i), 
where the constants a and b are defined by 



1 _ {m-aon){m-aon-ko) .^^ 
a =-r^ aon + m -\ _ ) > 0, 



:i.6) 



6 =1 + ^J^±^£ii^= > 0. 

\/{m + a^n — k^Y + Ak^aQU 

(ii) Uniqueness: Uniqueness holds m C (M+; (5^/2-1,^/2)1+1 ^ ^B'^ny) ifd^S. If 
d = 2, one should also suppose that rfiQ — rh, uq — n E j^^A+s g^j^^j^ ^^ ^ ^e jqj, q 
e G (0, 1), to get uniqueness m C(]R+; {B'^'^Y^^ x (B^Y). 

For the general data bounded away from the infinity and the vacuum, we have 
the following local well-posedness theory (one can refer to Theorem 3.1 for the corre- 
sponding statement in terms of new variables). 

Theorem 1.2 (Local well-posedness for general data). Let d^2,fl>0,2fi + dX^0, 
the constants m > and n ^ 0. Assume that itIq^ — fn^^ G iJ'^/2.'i/2+i^ f^^ _ f^ ^ 
j^d/2,d/2+i ^^^ ^^ ^ j^d/2~i,d/2 ^ j^ addition, sup^-g^d ?Tio(a;) < 00. Then there exists a 
positive time T such that the system (1.1) has a unique solution {rh, n, u) on [0, T] x M*^ 
and that {m~^ — fh~^,n — n,\i) belongs to F^ and satisfies s\iY)!tx)£\oT]y.m.<i''^iPi^) < ^^■ 

We also have the following continuation criterion for the local existence of the 
solution (see also Proposition 3.6). 

Theorem 1.3 (Continuation criterion). Under the hypotheses of Theorem 1.2, assume 
that the system (1.1) has a solution [rh, n, u) on [0, T) x W^ such that (m~^ —rh~^,n — 
n, u) belongs to F^, for all T' < T and satisfies 

m^'-m~\h-ne L°°([0, T); 5^/2,^/2+1)^ 

sup m{t,x)< 00, / ||Vu||ooC?t < 00. 

(t,a;)e[0,T)xRd Jo 

Then, there exists some T* > T such that {rh,n,u) may be continued on [0,T*] x R'' 
to a solution of (1.1) such that {rh~^ — rh~^,n — n,u) belongs to F^,. 

Remark 1.4. The results of the present paper are independent of the special structure 
(1.3) of the nonlinear pressure term P. Indeed, the similar results hold true as long 
as the term 'WP/fh can be decomposed into a linear term involving the modified mass 
and some nonlinear terms, similarly as (2.1) in the next section. 
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The rest of this paper is organized as follows. In Section 2, we investigate the 
global well-posedness of the Cauchy problem. We first reformulate the system through 
changing variables in order to obtain a priori estimates in the subsection 2.1. In the 
subsection 2.2, we are devoted to deriving a priori estimates for the transport equa- 
tion and the linear coupled parabolic-hyperbolic system with convection terms. The 
subsection 2.3 involves the proof of the existence and uniqueness of the solution. In 
Section 3, we prove the local well-posedness of the problem through some subsections 
similar to the global ones. An appendix is devoted to recalling some properties of the 
Littlewood-Paley decomposition and Besov spaces which we have used in this paper. 



2. Global well-posedness for small data 

2.1. Reformulation of the system. Let n ^ and m > (1 — sgnn)A;o, we introduce 
new variables n = aQ{h/m — n/m) and m = a{m — m) + bn, i.e. m = rh + [m — bn)/a 
in order to cancel the linear terms involving one modified mass from the momentum 
equation, where a and b are positive constants defined in (1.6). We also denote 
no = doino/fho — n/ffi) and rrtQ = a(mo — rh) + briQ throughout the sections for the 
global well-posedness theory. Then, we have 
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a/m 1 + ^— +nj m - /coj +AkQ\n + ^—\m. 



Cofh 

where the nonlinear term is 

XT/ X Vm-bVn { . _ _. . _^ ,,_ _nn\ 

H(m, n) := — — [aon + m)m + [am + b[aon + m))n 1 

a'^m^m \ / 

+ {K{m, n) - K(0, 0))| (n + — + l)mVn + (n + — + ifVm 

I m m 

aoh Vm — bVn 



(2.1) 



, , aon ^.\/m — b\/n , „ ^ 

koin +-^- 1) z + koVn \ 

m am ) 



N I / Cion , 'Wnim — bn) „ 

K(0, 0) Un + ^ + 1) ^ + rfin\/n 



m 



+ [n^ + 2n(-^ + 1) h -^niVm - bVn) 

m a am 

aon (m — bn)('Vm — bVn) 

-ko[— 1; 



m ' a^mfh 



Here, 

K{m,n) 



[(^ + m^) (n + ^ + 1) - A;o] ' + 4ko {n + ^) {m + ^) 



and K{0, 0) = Xj \/{m + aon — koY + Akoaon > 0. 



C. C. HAO AND H.-L. LI 



Therefore, with the new unknowns, we can rewrite the Cauchy problem of the 
system (1.1) as follows 



(2.2) 



where 



rij + u ■ Vn = 0, 

m^ + u ■ Vm + arhdivu = F{m, n, u), 

U( + u ■ Vu — /iAu — (yU + A)Vdivu + CoVm = G{m, n, u), 

{m,n,n)\t=o = (mo,no,Uo), 

/i = fl/rh, X = X/rh, F{m, n,u) = —{m — bn)dwu, 

G(m, n, u) = — CoH(m, n) z — (/iAu + (/i + A)Vdivu). 

am 

Note here that the first equation in (2.2) is a homogeneous transport equation, the 
estimates of n depend only on those of the velocity u. The second and the third ones 
in (2.2) consist of a coupled parabolic-hyperbolic system with the modified mass m 
and mixed velocity u involved. Thus, with the help of the decomposition (2.1), the 
original system is decoupled into a transport equation for the modified gas flow and 
a coupled system for the motion of the modified liquid fluid. 
We can get the following result for the reformulated system. 

Theorem 2.1. Let d ^ 2, n '^ 0, rh > [1 — sgnn)ko, fi > and 2fi + dX ^ 0, in 
addition, /j, + X > if d = 2. There exist two positive constants rj and Q such that if 
mo, uq E _B<^/2-1-'^/2 Qnd uo G B'^^'^~^ satisfying 

\\m,Q\\ ^d/2-l,d/2 + ||?T,0 ||_Brf/2-l.d/2 + ||Uo||5d/2-l ^ Tj , (2.3) 

then the following results hold: 

(i) Existence: The system (2.2) has a solution (m, n, u) in E'^^'^ which satisfies 

\\{m,n,u)\\Ed/2 < Q{\\mo\\Bd/2-i,d/2 + \\no\\Bd/2-i,d/2 + ||uo||Bd/2-i). 

It also belongs to the affine space 

(mi,no,Ui) + (Ci/2(M+;5'^/2"i))i+i x {C'/\R+; B''/^-'/')f, 

where (mi,Ui) is the solution of the linear system 

dtfTiL + arhdivuL = 0, 

dtUL - /iAuL - (/i + A)VdivuL + CoVmi = 0, (2.4) 

{mL,UL)\t=0 = ("^o,uo)- 

(ii) Uniqueness: Uniqueness holds in E'^^'^ if d ^ 3. If d = 2, one should also 
suppose that no,mQ G 5^'^+^ and Uq G B'^ for a e G (0, 1), to get uniqueness in E^. 

With the help of Theorem 2.1, we can prove Theorem 1.1 as follows. 

Proof of Theorem 1.1. From the conditions, we have uq G i?'^/2-i'^/2_ j^ addition, 
from ifiQ — ffiE iJ"'/2-i.a!/2^ ^g ^g^j^ derive mg = a(^o — fn)/m + bnQ G iJ'^/2-i.'i/2_ gj^ce 
(1.5) implies (2.3), the conclusion of Theorem 2.1 follows for {m, n, u). Changing back 
to the original variables (m,n,u), it leads to Theorem 1.1. By Lemma A. 3, it is easy 
to see that fh — m and fi — n also belong to C(]R"''; 5^/^"^ n B^^"^). 
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2.2. A priori estimates for linear system with convection terms. We first 
investigate some a priori estimates for tlie linear system with convection terms 
' rif + V ■ Vn = 0, 
m/ + V ■ V?Ti + amdivu = F, 

Ui + V ■ Vu - /iAu - (/i + A) Vdivu + CoVm = G, ^ ' ' 

^ (m,n,u)|i=o = (mo,no,Uo). 

We do not need to separate the velocity into the compressible and incompressible 
parts. In fact, we can prove the following proposition. 

Proposition 2.2. Lei a > 0, m > 0, s G (l-d/2, (i/2 + 1] andsi, S2 E {-d/2,d/2 + l] 
be constants. Assumes E L^{[0,T]; B''-^'^'^^) and denote V (t) = J^ ||v(r)||^d/2+i(ir. Let 
(m,n, u) be a solution of (2.5) on [0,T], then the following estimates hold: 

ll-nll- <-' oC'^(^)|Iti II (0 fi^l 

IPIIl°°([0,T];B''1'^2) ^ 6 |Po||B^1'''2, l^-Oj 



and 



I"^IIl°°([0,T];B''-i.'') + il^llL°°([0,T];B''-l) 



+ ||"^IUi([0,T];B»+i'=) + I|u||li([0,T];S'=+1) (2.7) 

<e (||mo||Bs-l,s + IIuoIIb"-! + ||-^||l1([0,T];B=-i.'') + ||G||ii([o_T];B»-i)) • 

Proof. Step 1: Estimates for the homogeneous transport equation. We derive the 
estimates for the first equation of (2.5) in Besov spaces. 

Applying the Littlewood-Paley operator A^ to the (2.5)-,^, it yields 

dtAkU + Afe(v ■ Vn) = 0, 

(2.8) 



Afen|j=o = Afcno. 
Taking the inner product of (2.8) with A^^n, we get for any Si, S2 G (— (i/2, 1 + d/2\ 

i-||Afen||H-(Afc(v-Vn),A,n)<7fc2^^'^^^'^^('=)||v||5,/2+i||n||Bn,»2||Afcn||2. 
It follows 

k J^ 

which implies the desired estimate (2.6) with the help of the Gronwall inequality. 

Step 2: Estimates for {m, u). Applying the Littlewood-Paley operator A^ to (2.5)2 
and (2.5)3, ^^ have 

(9jAfcm + Afc(v- Vm) +amdivAfcU = AfcF, 

(2 9) 
^iAfcU + Afe(v ■ Vu) - /iAAfcU - {fi + A)VdivAfcU + CoVA^m = A^G. ^ ' ^ 

Taking the inner product of (2.9)^ with Afcm and — AA^m, and (2.9)2 with A^u, 
we obtain 

- — \\Akm\\l + (Afe(v ■ Vm), A^m) + am(divAfcU, A^m) = (A^F, A^m), (2.10) 

--WVAkfnWl + (Afc(v ■ Vm), -AA^m) - am(divAfcU, AA^m) (2.11) 

= -(AfcF,AAfcm), 

^-||Afcu||2 + (Afc(v ■ Vu), Afcu) + /ill VA,,u||2 + (/i + A)||divA,,u||2 (2.12) 
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+ Co(VAfcm, Afcu) = (AfcG, A,.u). 
For the intersected term, we have 

-(AfcU, VAfcm) + (Afe(v ■ Vu), VAfem) - (Afc(v ■ Vm), divA^u) 

- am||divAfeu||2 + Co|| VAfcm,||2 + (2/i + A)(divAfcU, AA^m) '^^•^^) 

= - (AfcF, divAfcu) + (AfcG, VAfcm). 

Let 

a^ := ^11 Afcm||2 + || Afcu||^ + (^^ + ^)-^ ||y^^^||2 ^ 2A(AfcU, VAfcm). (2.14) 
am am 

For A = (/i + A) /(2am) > 0, there exist two positive constants ci and C2 such that 

cial ^ \\Akm\\l + II Afcu||2 + || VAfcm||2 ^ C2al (2.15) 

since we have, for M G {am/{2fi + A), 2afh/{fi + A)), that 

|2(AfcU, VAfcm)| ^ M||Afcu||2 + ||VAfcm||2/M. 
Combining (2.10)-(2.13), it yields, with the help of Lemma A. 4, that 

-— Ofc + /i||VAfcu||2 + (yu + A - amA) II divAfcull 2 + CoA||VA;,m||2 

- ^(A,F, A^m) + (A,G, A^u) - l^/i±^(AfcF, AA^m) - A^A^F, divA^u) 



am, am, 

+ A(AfcG, VAfcm) - ^(Afc(v ■ Vm), A^m) - (Afc(v ■ Vu), A^u) 

+ ^^^+_^^ (Afc(v ■ Vm), AAfcm) + A(Afe(v ■ Vu), VA^m) 
am, 

+ yl(VAfe(vVm),Afeu) 
<(|| AfeF||2 + II VAfeF||2 + II AfcG||2)(|| Afcm||2 + || Afcu||2 + || VAfcm||2) 
+ 7fc2"*^'^'"'"('=)||v||5,/2+i||m||B.-i,.||Afcm||2 
+ 7fc2-*^("-^)||v||Bd/2+i||u||B.-i||Afcu||2 
+ 7fc2-^('^^""('=)-i)||v||5./2+i||m||B.-M||VAfem||2 
+ 7fc||v||B./2+i(2-^(^-i)||VAfem||2||u||B.-i+2-'=(^^""(^^)-i)||m||B.-M||AfcU 

<(|| AfcF||2 + II VA,,F||2 + II AfcG||2 + 7fc2-'^'-^)||v||5./2+i(||m||B-M + HuH^.-i 
X (||Afcm||2 + ||Afcu||2 + ||VAfcm||2). 
Thus, it follows 
1 d 



^^«fc + Comin(22^1)a2 



which implies 

2'=(^-i)afc + Co / min(22^ l)2'=("-^)afc(r)c/r 
Jo 



G||^3-l + ||v||^d/2 + l I ||?Ti||5S-l,S + ||U||^3-1 



ttfc, 
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^0 



1^3-1, s + ||G||^3-i + ||v||^d/2+i ^^^2 '•'^ 'ak 



k 



Thus, by the Gronwall inequahty, we have 



I"^IIl°°([0,T];B''-1'=) + II^IIl°°([0,T];B''-1) 



Step 3: The smoothing effect for u. By (2.12), we have 

<||A,u||2(2^||Afcm||2 + ||AfcG||2 + 7fc2^'^'"^)||v||5./.+i||u||B.-i; 



It follows that 
d 



^ 2^'(^-i)|| Afeu||2 + C 5^ 2^'(^+i) II Afcul 

fc>0 fc^O 



<^2^-(^-i)[2^||A,m||2+||AfcG||2 + 7fe2"'^'^^)||u||5./.+i||ub.-i] 
<^2'''||Afem||2 + ||G||s.-i + ||v||5d/2+i||u||Bs-i, 

fc^O 

which implies, with the help of (2.16), that 

r^2^'(^+i)||A,u(r)||2 

<||uo||b-i+ / 5^ 2^=^11 Afcm(r)||2dr+ / ||G(r)b.-idr 

Jo ,.^n Jo 



+ ll"^l|Li([0,r];B»+l.'') + ||u||l1([0,T];B''+1'''-1) (2.16) 

<e (||?TIo||b''-i.'' + ||uo||b=-i + ||-^||l1([0,T];B=-i.'') + || G||l1([o,T];B''-1)) • 



+ sup ||u(r)||B.-i / ||v(r)||^d/2+idr 
Te[o,t] Jo 

<e^^W (||mob.-M + ||uo||B-iy" [||F(r)b.-i. + ||G(r)b.-i] dr") . 

Combining with (2.16), we get (2.7). D 

From the proof of Proposition 2.2, we immediately have 

Corollary 2.3. If a hounded operator B acts on the convection terms in (2.5), then 
the same estimates hold for the refined system 

ut + i3(v ■ Vn) = 0, 

rut + B{y ■ Vm) + amdivu = F, (2-17) 

Ui + S(v ■ Vu) - /iAu - (/i + A) Vdivu + CoVm = G. 
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2.3. Global existence and uniqueness of the solution. Step 1: Friedrich's ap- 
proximation. Let Lj be the set of L^ functions spectrally supported in the annulus 
Ci := {^ e W'- : 1/i ^ \^\ ^ i} endowed with the standard L^ topology. In or- 
der to construct the classical Friedrichs approximation, we first define the Friedrichs 
projectors (F^)^^^ by 

for any / G L^(M'^) where Ic^iO denotes the characteristic function on the annulus 
Ce. Then, we can define the following approximate system 

fn^ + F^(FV- VFV) = 0, 

mf + F^fF^u^ ■ VF^m^) + amdivF^u^ = F\ 

f 2 1 8") 

u^ + F^(FV ■ VFV) - /iAFV - (/i + A)VdivFV + CoVF^m^ = G^ 

where 

mi =F^mo, rii = FiUq, u^ = F^Uq, 

F^ =¥eF{¥em^, F^n^ F^u^), G^ = F^G(F£m^ F^n^ F^u^). 

It is easy to check that it is an ordinary differential equation in L^ x L^ x {LjY 
for every £ e N. By the usual Cauchy-Lipschitz theorem, there is a strictly positive 
maximal time T/ such that a unique solution (m^,n^,u^) exists in [0,T/) which is 
continuous in time with value in Lj x Lj x {LjY, i.e. {m^,n^,u^) G C([0,T/);L^ x 
Lj X [LlY). As F^ = F^, we see that F£(m^, ra^, u^) is also a solution, so the uniqueness 
implies that Ff(?n^,?T,^, u^) = {m^,n^,u^). Thus, this system can be rewritten as the 
following system 

mf + F^fu^ ■ Vm^) + amdivu^ = Ff, 

('2 19) 

^ (m^ra^u^)|i=o = {mi,ni,ui), 
where 

F/ = ¥i,F{m\ n\ u^), and G^^ = ¥eG{m\ n\ u^). 
Step 2: Uniform estimates. Denote 

Eo = ||"^o||Bd/2-l,d/2 + ||?^o||Bd/2 + ||Uo||Bd/2-l, 

and 

T, := sup{T G [0,r;) : \\{m\n\u')\\^,j2 ^ ACE^], 

where C corresponds to the constant in Proposition 2.2 and A > max(2, 1/C) is a 
constant. Thus, by the continuity, we have Ti > 0. 

Let Mq be the continuity modulus of the embedding relation 5*^/^ (W^) ^-> L^ 
We make the assumption 

2(1 + b)ACMoEo ^ am. 

Then, it implies 

11^ ||L°°([0,T]xRrf) ^^o||"^ IIl°°([0,T];_B'*/2) ^ ^o||"^ llL°°([0,T];_B<^/2-i.d/2) 



' oo/ 
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afh 



^ACMoEo ^ 



2(1 + 6)' 
Similarly, we have 



Then, 



\ i] am 
\n |L-([o,T]xRd) ^ 2(1 + 6)' 


^ , f fs m — bn 
m,[m , n ) = m -\ G 

(Ji 


m Sm 

2' 2 



By Proposition 2.2, we have 



(m^n^u^)||5^ 



From Lemmas A.l and A. 3, we get 

II T?^\\ 
lKlllLi([0,T];Bd/2-i.d/2) 

^\\m dlVU ||^i([o,T];B<'/2-i>'*/2) + IP dlVU ||^i([o,T];Bd/2-i.d/2) 

^(11"^ IIl'«([0,T];B'*/2-i.'*/2) + \\n ||ioo([o,2.].Bd/2-i,d/2))||divU ||Li([0,T];Bd/2) 

^{\\m ||ioo([o,T];B'*/2-i.'*/2) + If llL^([0,T];Bd/2-i.rf/2)j||U IIl1([0,T];B'*/2+1) 



<l 



772 ,n , U 



|2 



Thus, 



)"- 5 " ;ii^d/2 



-(/iAu^ + (/i + A)Vdivu^)||^i([o_T];Bd/2-i) 



^^ll£ll ll£ll Njl^ll 

i^llF llL°o([0,T];Bd/2) + ll""^ llL°°([0,T];Bd/2)j||U llLi([0,T];Bd/2+l) 

<||(m^n^uO|||,. 
Similarly, we can get 

||H(m ,n )||Li([0,r];_Bd/2-i) <{\\n ||L°°([0,T];Bd/2) + 11"^ llL°°([0,T];Bd/2)) 

<||(m^n^uO|||./.. 
Hence 

^Ce^'^^«(l + CA2C''Eo)^o. 
Thus, we can choose Eq so small that 

1 + CA^C^Eo ^ -^—, e^'^^« ^ ^^ and 2(1 + 6)A(7MoEo ^ am, (2.20) 

which yields ||(m^n^ uO||^d/2 ^ ^ACE^ for any T < T^. 

We claim that T^ = T/. Indeed, if Ti < T/, we have seen that ||(m^,n^, u^)||^d/2 ^ 

^^^AC'E'o. So by the continuity, for a sufficiently small constant s > 0, we can obtain 
II {m,^, n^, u^) II d/2 ^ ACEq which contradicts with the definition of Ti. 

(T + a) 
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Now, we show the approximate solution is a global one, i.e. T/ = oo. We assume 
T/ < cxD, then we have shown \\{m^,n^,u^)\\^d/2 ^ ACEq. As 

m' G L~([0,T;);5'^/2-i,d/2)^^£ ^ L~([0,T;); ^-^/^^ and u^ G L~([0,T;); E'^/^-i)^ 

it implies that 

||(mV^u^)||L-([o,^;);L|) < oo- 
Thus, we may extend the solution continuously beyond the time T/ by the Cauchy- 
Lipschitz theorem. This contradicts the definition of T/. Therefore, the solution 
(m^, n^, u^)i^n exists global in time. 

Step 3: Time derivatives. For convenience, we split the approximate solution 
(m^,n^,u^) into a solution of the linear system with initial data (m£,n^, u^), and 
the discrepancy to that linear solution. More precisely, we denote by (r7i£,uf^) the 
solution to the linear system 

ml + amdivu^ = 0, 

u^ - /iAu^ - (/u + A)Vdivu^ + CoVm^ = 0, (2.21) 

and (mfp, n^, u^) = (m^ — m^, n^ — rii, u^ — u£). 
It is clear that the definition of {mi,ni,U£) implies 

me ^ mo in E^/^-^'^/^^ n, ^ rio in B^^'-'^^/' , u, ^ uq in 5^^/^"^ 

From Corollary 2.3, we have 

{m[,ne,u1) -^ {mL,no,UL) in £"^/^, 

where m,L and u^ satisfy the linear system (2.4). 

Now, we derive the uniform boundedness of the time derivatives of the discrepancy 

(m^,n^,u^). 

Lemma 2.4. ((?72^,nf), u^))^gN is uniformly bounded in 

Proof. Since 

9y^ = -F,(u^VnO, 

we have ^jn^ G L2(M+; E^/^-i^ gj^^^g ^£ ^ L°°(M+; B"^/^) and u^ G L2(M+; fi'^/2) with 
the help of the interpolation theorem. 
From the equation 

dtm^u = — F^(u^ ■ Vm^) — amdivu^ + amdivu^ — ¥i{{m^ — 6n^)divu^), 

it follows that dtm'^^ G L'^{R+; B'^l^~^). 
Recall that 

a^u^ = - F^(u^ ■ Vu^) + /iAu^ + /iAu^ + (Ai + A)Vdivu^ - (a* + A)Vdivul 

-CoVm^-CoVml + G^i, 

we can obtain ^tuf, G (L°° + L®/^ + L®/^)(M+; E^'/^"^/^) through easy but tedious 
computations with the help of Lemmas A.l, A. 3 and A. 5. 

Applying the Morrey embedding relation ly^'^(M) C C^^^/^(M) to the time variable 
for 1 < p ^ oo, we obtain the desired result. D 
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Step 4: Compactness and convergence. The proof of the existence of a solution is 
now standard. Indeed, we can use Arzela-Ascoh theorem to get strong convergence 
of the approximate solutions. We need to localize the spatial space in order to utilize 
some compactness results of local Besov spaces (see [1, Chapter 2]). Let (Xp)peN be a 
sequence of ^(M'^) cut-off functions supported in the ball 5(0, p+1) of M*^ and equal to 
1 in a neighborhood of B{0,p). In view of Lemma 2.4 and uniform estimates obtained 
in Step 2, we see that ((Xp^z?? Xp^d? Xpuf)))^6N is bounded in E'^^^ and uniformly 
equi-continuous in 

C([0,T];(5^/2-i)i+ix(5'^/2-5/4)'i) 

for any p eN and T > 0. Moreover, the mapping / H- Xpf is compact from 5'^/2-i>d/2 

into 5^^/2-1 ^^J f^Qj^ ^d/2-l ij^^Q 5^/2-5/4 _ 

Applying the Arzela-Ascoli theorem to the family {{Xp^^Dy Xp'^^d^ Xpui)))^eN on the 
time interval [0,p], then we use the Cantor diagonal process. This finally provides 
us with a distribution {vfiD^nD^'^D) continuous in time with values in (i^'^/^-i^i+i x 
(^5<^/2-5/4^d and a subsequence (which we still denote by the same notation) such that 
we have for all p G N 

(Xp"^l», Xp^D^ XpU^) -> {Xp'^D, Xp^D, Xp^d), as £ -> cx) 
inC([0,p]; (5'^/2-i)i+i x (5^/2-5/4^^^^ r^^^^ obviously implies that (mfj,?!^^, uf^) tends 
to {mD,nD,UD) in ^'{R+ x M^). 

Coming back to the uniform estimates and Lemma 2.4, we further obtain that 
{mD,nD,UD) belongs to E'^/'^ and to (Ci/2(^+. 5^/2-1)^1+1 ^ (ci/8(k+. 5d/2-5/4))d_ 

The convergence results stemming from this last result and the interpolation argument 
enable us to pass to the limit in ^'{M.^ x R'') in the system (2.18) and to prove that 
(m, n, u) := (m^,, ul, ul) + {ma, no, uo) is indeed a solution of (2.2) with the initial 
data. Since it is just a matter of doing tedious verifications, we omit the details. 

Step 5: Continuities in time. The continuity of u is straightforward. Indeed, from 
the third equation of (2.2), we have Uf G (L^ + L^)(]R+; (i?''/^"^)'^) which implies 
u G C(M"'"; (i^'^/2-i)'') in view of the Morrey embedding and the embedding relation 
ly^'^(M) C C(]R). Consequently, the continuity of n in time is obtained from (2.6). 
For m, it is easily to see that m^ G L^(M+; 5'^/^"-'^) fl L^(M+;i?'^/^) from the second 
equation of (2.2) which yields m G (^(M"*^; 5'^^^"^''^/^) by the embeddings mentioned 
above. 

Step 6: Uniqueness. Next, we prove the uniqueness of solutions. Let (mi,ni,Ui) 
and (7712, ?T-2, U2) be two solutions of (2.2) in Ej! with the same initial data. Denote 
((5m, (5n, 5\i) = (m2 — mi, 77-2 — ni, U2 — Ui). Then they satisfy 

dtSn + U2 • V6n = —6u ■ Vrii, 

dtSm + U2 ■ V6ni + amdiv^u = —6u ■ Vrui + 5F, 

(2 22) 
dt6u + U2 ■ V6u - /iA(5u - (/i + A) Vdiv^u + CoVSm = -6u ■ Vui + 6G, ^ ' ' 

[(5m,5n,5u)|i=o = (0,0,0), 

where bF = F(m2, n2, U2) — -F(mi, rii, ui) and 6G = G(m2, n2, U2) — G(mi, ni, ui). 
We first consider the case d ^ 3. Similar to the derivation of (2.6), we can get for 
t G [0, T] with the help of (2.6) 

C C^ II II fl 

P'^llL°°([0,r];B'*/2-2.d/2-l) ^e ll"2|l3d/2+l ^||5u||^l([o^y].3d/2)||ni||icx>([o,T];Bd/2-l,d/2) 

<e^^""^"^'"^^"^'*/^+^'ll^o||B^/2-i,^/2||5u|Li([o,T];B./2). (2.23) 



2-1, 
T 
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By Lemma A.l and A. 3, we have 

<e^fo'^\\^^h^/^+^''-[\\6u-Vm4LH[o,T];B^/^-^^'^/^-^) + PF(r)|Li([o,T];i?<*/2-2../2-i) 

+ \\6u ■ VUi||^i(jo^^].^d/2-2) + 11(56(7") ||ii([o,T];Bd/2-2)] 
<e^^"'""'"^'/'+''^"(pu||il([0,T];B<^/2)||mi||^^(P^^].5,/2-l.d/2) 

+ ||(5"^,5n)||^«=Qo^j,].^d/2-2,d/2-i)||u2||ii([o,T];Brf/2+i) 
+ ll("^l)'^l)llL°°([0,T];S''/2-i^d/2)||5u||^i([0,T];B'«/2) 

+ ||5u|lz°°([0,T];Brf/2-2)ll"l|lLi([0,T];Bd/2+i) + PG(t) ||l1([0,T];B'^/2-2) 
<e^il"2ll,iao,r,../2+i,(^(l + ||K,n2)|k.([0,T];B^/2-M/.)) 

X ll("^i,^i)llL-([o,T];B''/2-i.d/2) + ^(T)j||((5m,(5n,(5u)||^d/ 

where hmsupj._^Q+ ^(T) = 0. Thus, 
\\{6m,6n,6u)\\j^d/2-i 
^Ce^\\i^-''^)h.ao^T,s^^..., ((1 + ||(m2,n2)||zo.([o,T];i../2-i../2)) 

X ||(mi,?7,i)||^oc([o,T];i?d/2~i.d/2) +^0 + ^(r)j||(^"^,^'^,MIU^/2-l- 

We take Eq small enough such that it satisfies the condition 2C(1 + CACEq)ACEq + 
Eq < 1/4 and (2.20), and choose T > so small that C*||(ui, U2)||^i([Q2-].5d/2+i) ^ ln2 
and Z(T) < 1/2, then it follows \\{6m,6n,6u)\\^d/2-i = 0. Hence, (mi, rii, ui)(t) = 

{m2,n2,U2){t) on [0,T]. By a standard argument (e.g. [5]), we can conclude that 
(mi,ni,ui)(t) = (m2,n2,U2)(t) on M+. 

For the case d = 2, we have to raise the regularity of the spaces. Thus, we also 
suppose that mo, no G 5'^'^+^ and uq G -B*^ for a e G (0, 1). By the same process, we 
can prove the existence of solution (m, n, u) in the space E'^^'^ provided the norms of 
initial data is sufficiently small. Then, in the same way as in the case d ^ 3, we may 
prove the uniqueness of solutions in the space E"^ (of course, holds in E^). We omit 
the details. 

3. Local well-posedness for large data 

3.1. Reformulation of the system. We change variables to p = ffi{fh~^ — ffi~^) 
and g = n — n. Then we can reformulate the system (1.1)-(1.2) as 

' Pi + u ■ Vp = (p + l)divu, 

gt + M-Vg = -{g + fi)divu, 

Ui + u ■ Vu - (1 + p)(pAu + (p + A)Vdivu) + Q(p, g) = 0, 
, (p,5',u)|i=o = (po,fi'o,uo), 

where po = m(mQ'"'^ — m"^), go = Uq — n and 

Q(p, g) :=m-i(l + p)VP(m/(l + p),g + n) 



(3.1) 
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pVp _ , Oo 



-Vp + ^{p + l)Vg + Bip,g) - — ^Vp + :^^^-:F Vp 



P + 1 TTl 



m ko — aon, 



{p + iy p + i 



+ — -^ + —{9 + gp)^g + = V^ + pVg 

p + 1 m m m . 



with 



r Tfi 9 

B{p,g):= ( — —r + ao{g + n)-ko) + 4:koao{g + n) 
I P+l 



-1/2 



'P 

We now state the result for the local theory for general data bounded away from 
the vacuum as follows. 

Theorem 3.1. Let d ^ 2, p > 0, 2p + dX ^ 0, the constants rh > and n ^ 
0. Assume that po G 5'^/2,d/2+i^ ^^ ^ ^d/2,d/2+i ^^^ ^^ ^ ^d/2-i,d/2_ j^ addition, 

inf Pq{x) > —1. Then there exists a positive time T such that the system (3.1) has a 



unique solution (p, (7,u) on [0, T] x R'^ which belongs to 

(C([0,T];5^/2,d/2+i))i+i ^ (C([o,T];5''/2-i.''/2) nLi([0,T];5^/2-fi,d/2-f2))d^ 

and satisfies inf p(t,x) > — 1. 

(t,x)e[0,T]xR'* 

3.2. A priori Estimates. Now, let us recall some estimates for the following para- 
bolic system which is obtained by linearizing the momentum equation 

ut + V ■ Vu + u ■ Vw - b{t, x)(pAu + (p + A)Vdivu) = /, 

which had been studied in [1, 8]. Precisely, we have the following lemma (cf. [1, 
Proposition 10.12]). 

Lemma 3.2. Let a G (0, 1], s G (— (i/2, d/2], u = min(p, A + 2p) and v = p+ |p + A|. 
Assume that 6 = 1 + p with p G L~([0, T]; E'^/^+a) ^^^ ^/^^^ 

6* := inf hit.x) > 0. 

(t,x)G[0,T]xR<' 

There exist a universal constant k, and a constant C depending only on d, a and s, 
such that for all t G [0, T], 

^ (IIuoIIb'' + ||/||Li([0,t];B»)) 

X exp (c j (||v||5./.+, + ||w||5./.+, -f- (M'''/"^'/llp|lB.V)rfr') . 

// V and w depend linearly on u, then the above inequality is true for all s G (0, d/2 + 
a] , and the argument of the exponential term may be replaced with 

C j' (llVulU + (M'"'/"^'/llp|l|'.V) dr. 

For the mass equations, we only need to study the following equation with two 
constants 9 &M. and /3 > 

h^ + w■Vh = e{h + /3)divv, 

h\t=o = ho. ^^-^^ 
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Proposition 3.3. Let s e {-d/2,d/2 + 1], T > 0, 9 e R and f3 ^ be constants. 
Assume that ho G 5''/^, v e L^([0,T); 5'^/^+^) and a satisfies (3.2). There exists a 
constant C depending only on d such that for all t G [0,T], we have 



Z°°([0,i];B'«/2) 



^gC(l+2|e|)/„'||v||^,/,+,d. (||, II 



/3 



/3 



l + 2\e\J l + 2\9\ 



(3.3) 



and 



Z°°([0,t];Bs) 



<eC(i+[e|)/J||v||^,/2+idr/||^^ 



I'OWb^ 



+ C'|^|[e^(i+2|^l)/olMI,./.+i'^-(||/,Q||^,/,+ 



/3 



1 + 21^17 1 + 21^1 



|v||ss+i(ir ]. 



(3.4) 



Proof. Applying the operator A^ to (3.2) yields 

dtAkh + Afc(v • V/i) = eAkiih + /3)divv). 

Taking the L^ inner product with A^h, we get, with the help of Lemmas A. 3 and A. 4, 
that 

1 C? II A , Il9 

- — \Akh\l 
2dt" "^ 

= - (Afc(v ■ Va), Afc/i) + ^(Afc((/i + /3)divv), A^/i) 

<7fc2-^-^||v||5./2+i||/i||BH|A,/i||2 + |^|7fc2-''^(||Mivv||B. + /3||divv||BOI|Afc/i||2 
<7fc2"^-^((l + m\h\\BsMs.r2+^ + (||/i||b^/2 + /3)\e\MBs+i)\\Akhh. 
Eliminating the factor ||Afc/i||2 from both sides and integrating in the time, we have 
\\Akh\\2 ^WAkhoh 

+ C7fc f 2-''{{l + \e\)\\h\\BsMB./2+^ + {\\h\\Ba/2+/3)\e\MBs+^)dr. 

Jo 
It follows, for any A; G Z and any t G [0, T], that 

2'''\\Akhh^2'''\\Akho\\2 

ft (3.5) 

+ C-fk / ((l + |^|)||/i||BH|v||Bd/2+i + (||/i||Bd/2+/3)|^|||v||B^+i)cir. 
Jo 

Summing up on A; G Z, it yields 

II^IIl°°([0,T];B») ^II^oIIb'' 

+ / C[{l + \e\)\\v\\Bd/2+i\\h\\Bs + {\\h\\Bd/2+(3)\e\\\v\\Bs+i]dT. 

Jo 
By the Gronwall inequality, we have (3.3) for s = d/2 and then for any s G 

(-rf/2, d/2 + 1] 



\L°°{[0,t];B=) 

^e^(i+l^l)/ollv|l,./2+i^-(^||/,^||^^ 



c\e\ 



,ca«.,/.-,M,_*(||,„||^,,,,_jL_),_M^^ 



v\\pis+idr 
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This completes the proofs. D 

In general, for the transport equation 

jht + ^r■Vh = f, 
I h{0) = ho, 
we can get, in a similar way with Proposition 3.3, that 
Proposition 3.4. Let Si, Sa G (-rf/2, d/2 + 1] andT >0. Then it holds for t E [0, T] 

\\h\\L^im;Bn,s,) ^ e^^olMU./2+irf- (\\ho\\Bn,s,+ f WfWBn-.s^dr 

3.3. Existence of local solution. Step 1: The Friedrich's approximation. For con- 
venience, we introduce the solution uis to the linear system 

9tUis - yuAuis - (/i + A) Vdivuis = 0, uis(O) = uq. (3.6) 

Denote uf^ := F^uis and u^ := u^ — uf^. Then we can construct the following approxi- 
mation {p^,g^,u^) satisfying 

f pf + F,(u^ ■ V/) = F,((/ + l)divuO, 

g', + F,(u^ ■ Vg') = -FeUg' + n)dwu'), 

dtu' + F,(uf3 ■ Wu') + F,(ii^ ■ VuO - F,[(l + p')ifiAu' + (/i + A)VdivuO] (3.7) 
= F4/(/iAuf, + (/i + A) Vdivufj] - F,(uf3 ■ Vufj - F,Q(/, /), 
, {p^,g^,u%=o = {pl9o,0), 

where p^ := F^po, Qo '■= ^i9o and u^ := uf^ -h u^ 

Note that if H-po is bounded away from zero, then so is 1-l-F^po ^ot sufficiently large 
i. It is easy to check that (3.7) is an ordinary differential equation in Lj x L'j x {LjY 
for every £ G N. By the usual Cauchy-Lipschitz theorem, there is a strictly positive 
maximal time T/ such that a unique solution {p^,g^,u^) exists in [0,T/) which is 
continuous in time with value in LjxLjx {LjY, i.e. (p^, g^, u^) G C([0, T/); LjxLjx 
{LjY), and 1 -|- p^ is bounded away from zero. 

Step 2: Lower bound for lifespan and uniform estimates. We introduce the following 
notations 

Mo :=||po||Bd/2.d/2 + Q, M (t) := \\p ||^oo([0^i];Bd/2,d/2 + a), 

^0 :=||5'o||Bd/2,d/2+Q, A^ (t) := \\g ||^oo([o^i];Bd/2,d/2+a), 

|Uo||Bd/2-l.d/2-l+Q, f^ls(^) '■= l|Uls||Ll([0,t];Bd/2+l,d/2+l+a-), 



I" llL°°([0,t];B''/2-l^'*/2-l + «) +0*Z/||U ||^l([o,t];_Bd/2 + l,d/2+l + a). 



Uo 
U\t) 
In view of Lemma 3.2, we take v = w = / = and p = there to get 

II ^ II < TT 

ll^lsllL°°([0,t];B''/2-l.d/2-i+<^) ri t^O- 

Note that for all /c G Z, we have 

||AfcPo||2 ^ IIAfcPolh, WPoWb- ^ IIpoIIb'' 

and similar properties for gf^ because of the boundedness of the operators A^. 
From (3.3) and (3.4), we get 
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\\P llL°°([0,t];B''/2+°=) 

+ c 



^ e 
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e3C/o'||u^lU.,...^||^^||^^^^+l) + 2jy^ 



|U \\j^d/2 + l + adT j , 



and 



II i>u ■^r r* ih/ii , riT /ii II '^A '^ 



(3.9) 
(3.10) 



+ C 



q^rt|i,/ii , j^/,, ,, n\ 2n1 /"^ 



\\9 llL°°([0,t];B''/2+") 



|U ll^d/2+l + adT I . 



(3.11) 



- J ^ u 

Let b{t, x) = 1 + p^. From Step 1, we know 6* > 0. Thus, by Lemma 3.2, A.l and 
A. 3, we have 

1^ llL«=([0,t];B<*/2-i) + K,b^K\\U llLi([0,t];B''/2+i) 



< 



ks)\\Bd/ 



+ llufs ■ VufJlBd/2-1 + ||(5(p^/)||Bd/2-l)rfr 

X exp fcj^ (||u1^./2+i + ||ufJ|^./2+i + (M'"'/"^'/ll/ll|'.V)c?r') , 
and similarly 



(3.12) 



^ y (||p'(/iAuf, + (/i + A)Vdivu[J||5,/2-i+. 

I 11..^ VT..^ II , 11^/ „£ „^\ 



-1+0 )dT 



I /L , Al-2/a-2/a|| ^^||2/a \ 7 
1 + [b^K) ' V ' \\P \\Bd/2+cjd'T^ 



(3.13) 



X exp ( C / M|u^||^d/2+i + IJufJI^d/ 
By Lemma A. 3, we get, for all a G {0, a}, that 

■■ l \ l W ^ W t\\ II A ^ II ^11 ^11 II f. II 

^U1sIIb<*/2-1+<^ ri IIP |lBd/2||^Uig||5d/2-l+<7 ^ \\p || 5d/2 || Ujg || 5d/2+l+<T , 



I „^Y7 J- ...^ II < 11/1^11 11,1^11 

IP vaivU[g||^d/2-i+<T ^ IIP ||^d/2||U[g||^d/2+i+<T, 

^ Y7,,^ II < 11,/ II IIY7,,^ II < 11,,^ II 11,,^ II 

1= ■ vuiji Rd/2-i+<T ;:^ ||U[g||^d/2|| vU[g||^d/2-i+<T ;:^ ||U[g||^d/2||U[g||^d/2+o 



11,/ Y7,,^ II < 

II "is ■ vuij|^d/2-i+<T ;i 
Recall that 

-' "'' =4^ - V/ + ^(/ + 1) V,^ + B{p', g') \ - ^T^V/ 



Qip',9' 



+ 



ko-apn , oofi^^V// a| , ^ ^ ^ ao(/i;o + m + apn) 
p^ + 1 p^ + 1 rh m 



p^ + 1 p^ + 1 m 



— ^^p'-Wg 

m 



where 



B{p',g'):= 



m 



p' + l 



\2 -| -1/2 

- ao(5''^ + n) - fcoj + 4:koao{g'^ + n) 
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Similarly, for the third term of Q 

IK/ + 1)V/||b./2-i+. < (II/IIb^/2 + l)||V/||5./2^i+. < (||/b./2 + 1)\W\\b^/^+^- 
By Lemmas A. 3 and A.l, we get for the first two terms of Q 



p' + l 



Vp^ 



< 



5d/2-l + (T 



p 



p' + l 



Bd/2 



1 II v/ 



|5d/2-l + CT 



<(||p^||Bd/2+l)||plBd/2+. 



and 



fi(/,/ 



V/ 



(/ + 1) 



< 



[i?(/,/)- 5(0,0)] 



^d/2-l + (T 

v/ 



(/ + 1)^ 

<(||E(/,/)-i?(0,0)||^d/2 + l) 



+ 5(0,0) 



^d/2-l + CT 



V/ 



(/ + 1)^ 



^d/2-l + CT 



p'p' + 2p' 



(/ + 1)^ 

;$ (Up llBd/2 + US' 11^^/2 + 1) ||p ||Bd/2+<T. 

Similarly, 



^d/2 



1 ||V/||sd/2-l + „ 






Vp' 



p' + l 



Qd/2-l + a 



^ / \\ i\\ ii£ii i\^ii-^ii 

^ (Up llBd/2 + US' llBd/2 + 1) \\p \\Bd/2 



p' + i 



J^d/2-l + a 



~ (IIP llBd/2 + 115' llBd/2 + i) \\g 



_Bd/2||P ||5d/2+<T, 
2 



II-d(,p ,g )[i + p )g yg ||Bd/2-i+a ^ (||p ||^d/2 + \\g ||^d/2 + i) \\g ||Bd/2||5' ||Bd/2+<T, 

ll-DlP ,9 )Vg llBd/2-l+.T ^ (IIP llBd/2 + 115 llBd/2 + ij 115 llBd/2+<T, 

II-d(p ,g)p\/g ||5d/2-i+a ^ (Up ||Bd/2 + \\g \\Bd/2 + i) \\p ||Bd/2||5' ||Bd/2+a. 
Thus, we get 

II Q(P^50 II i?d/2-l+. < (||plBd/2 + ||5^||Bd/2 + 1)^ (||p1b<^/2+<. + 115^5^/2+0 • 

Therefore, from (3.8)-(3.13), we conclude that 
M\T) <eCmT)+uHT)/ib,u))^^^ ^ 1) 

+ e^(^^(^)+'^^(^)/(''*^))(Aro + f )(f/,l(T) + f/^(T)/(M) - |, 
U\T) < {{uM\T) + Uo)Ui{T) + {M\T) + N\T) + lf{M\T) + N\T))T) 

Now, if we take T so small that 

exp {CUi{T)) ^ 72, exp {^^^ 1 ^ v^, 
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and 

exp (C(M'"'^"'^'^"(^^'(^))'^"r) ^ 2, 
then we have 

v;^ 0^3, v;^ °^3' (3.14) 

U'{T) ^ C ((Mo + iVo + 1)'(T + i^UiiT)) + UoKiT)) . 

Noticing that {¥ip^,¥£g^,¥^u^) = {p^,g^,u^) by the construction of the approxi- 
mated system. Thus, we have 

dtil + pY' + ^ii^' ■ V(l + pY') ± ^^(1 + p')^divuO = 0. 
It follows, by noticing that \dt\f\\ = \dtf\, that 

11(1 + /)^^WI|oo ^11(1 + P^)^^||oo 

+ Allu^' ■ V(l + /)^ioo + 11(1 + /)^divu^lloo]rfr, 
which yields, by the Gronwall inequality, that 

||(l + pO±i(t)|U^e/oll<i'™^ll-'^-('||(l + p^J±i|U + |*||u^|U||V/|Urf^^ 

^efo\\'^-^^-'\\-^^(^\(l+pX'\\oo + cl\\\ui\\^.;. + \\u^^^ 
^e^o lldWlloorfr A(i + p^)±i||^ + CT{Uo + U\T))M\T) 



where we have to choose a = 1 in the previous estimates. Hence, if we assume that 
there exist two positive constants 6* and h* such that 

then we can take T small enough such that 

lldivu^lloorfr ^ In 2, and CT{Uo + U\T))M\T) ^ 1, 



^ 1 + p^ ^2(6* + l). (3.15) 



'0 
and so 



2(1 + 6J 



Now, by means of a bootstrap argument, we can get that there exist two constants 77 
and C depending only on d such that if 



(Mo + iVo + 1)\T + uUiiT)) + UoUiiT) ^ rjKu 



4,^, -..,^.. . .....^. . , (3-16) 



then we have (3.14) and (3.15). 

Therefore, T/ may be bounded from below by any time T satisfying (3.16), and the 
inequalities (3.14) and (3.15) are satisfied by {p^,g^,u^). In particular, {p^,g^,u^)i^fq 
is bounded in F^. 

Step 3: Time derivatives. In order to pass to the limit in the approximated system, 
we first give the following lemma. 
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Lemma 3.5. Let p^ := p^ — F^po? 9^ '■= 9^ — '^1-90- Then the sequences {p^)een o,nd 
{g^)een are uniformly bounded in 

C([0, T]; 5'^/2.'i/2+i) n Ci/='([0, T]; 5'^/2-i,rf/2)^ 
and the sequence (u^)^^^ is uniformly bounded in 

(C([0, Tl; 5^/2-Ld/2\ p C^/^HQ^ T]; 5'^/2-l.'i/2 _,_ ^d/2-3/2,d/2-l/2NN<i_ 

Proof. From dtp^ = — F^(u^ ■ Vp^) + ^e{{p^ + l)divu^), we have 

1 1 o ~£ 1 1 ^ 1 1 A ~^ 1 1 

ll^tP llL2([o,T];B''/2-l.d/2) ^ ll'^tP ilL2([0,r];_B''/2-l.d/2) 

^11" ■ '^P llL2([0,T];Bd/2-i,d/2) + ||(p + l)divU ||£2([o^r].Bd/2-l.d/2) 

<ll ^11 II ^11 

~ll" llL2([0,T];Bd/2)llP llL°°([0,r];Bd/2,d/2+i) 

+ (Up llL°°([0,r];_Brf/2) + 1)I|U IIZ2([0,T];B'*/2,d/2+i). 

Since (u^)£gN is uniformly bounded in 

it is also bounded in L^([0, T]; 5'^/^''^/^+^) by Lemma A. 5. Recall that {p^)e<zn is uni- 
formly bounded in Z°°([0,T]; 5'^/2.^/2+i^^ then ((9(P^)^gN is uniformly bounded in 

L2([0,T];5'^/2-i,d/2)^ 

and so (p^)£6N is uniformly bounded in 

C^/'([0,T];5'^/2-i,d/2^ g^^^ i^ C([0,T];5'^/2,<i/2-i^)_ 

Similarly, we have the same arguments for g^. 
Recall that 

dtu^ = - ¥^{u^ ■ VuO + F4(l + pO(pAG^ + {fi + A)VdivGO] 

+ ¥e[p\fiAui + (/i + A) Vdivufj] - F,Q(/, /). 

Since 

||u • Vu ||^2([o,T];Brf/2-i.d/2) < ||U ||/,oo([o,T];Bd/2-l.d/2)||u ||L2([o^r].Bd/2+i) , 
11(1 + p^)(yUAu^ + (P + A) Vdivii^) ||^4/3([o,T];iJd/2-3/2.d/2-l/2) 
^(1+ Up llL°°([0,T];_Bd/2))||u ||^4/3([o,T];Bd/2+i/2,d/2+3/2), 

||p^(pAuf, + (p + A)VdiVufj||i4/3([o,T];i?<*/2-3/2.d/2-l/2) 

^ II £|| II ^ II 

^IIP llL°°([0,T];Brf/2)l|UlsllL4/3([0,T];Sd/2+i/2,d/2+3/2), 

l|Q(P 'fi' )llz°°([0,T];B<*/2-i.'i/2) 

^ (Up llz°°([0,T];Brf/2 + US' llz°°([0,T];Brf/2) + ^ ) M|P llz°°([0,T];Bd/2,d/2+i) 

II ^11 ^ 

+ 115^ IIZ°°([0,T];B'«/2.d/2+i) ) 7 

by Lemma A. 5, we can conclude that (u^)^^^ is uniformly bounded in 

C'/\[0,T]; 5'^/2-i.'^/2 + 5<i/2-3Ad/2-i/2) ^^^ .^ C([0,T]; 5'^/2-i,d/2)_ 

This completes the proof of the lemma. D 
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Step 4- Compactness and convergence. The proof is based on the Arzela-Ascoh the- 
orem and compact embeddings for Besov spaces. Since it is similar to the arguments 
for global well-posedness, we only give the outlines of the proof. 

From Lemma 3.5, {p^)eeN is uniformly bounded in the space 

Z°°([0,T];5^/2,d/2+i) 

and equicontinuous on [0,T] with values in 5'^/2-i,'^/2_ gj^ce the embedding 

^d/2-l,d/2 ^ Qd/2-l 

is (locally) compact, and (po)£gn tends to po in i?'^/2.''/2+i^ ^g conclude that {p^)ii^n 
tends (up to an extraction) to some distribution p. Given that (p^)£gn is uniformly 
bounded in L~([0,T]; 5*^/2 j^ we actually have 

peL^{[Q,T]-B'^'^). 

The same arguments are valid for the sequence {g^)i^n- 

From the definition of (ufj^g^, it is clear that (ufg)£gN tends to the solution uis of 
(3.6) in L~([0,t];5'^/2-M/2) p^iQO,T] ; 5^^/2+^/2+2)^ 

Since (u^)^gN is uniformly bounded in L°°([0, T]; i?'^'^"^''^'^) and equicontinuous 
on [0,T] with values in 5'i/2-i,'^/2 _^ 5^/2-3/2,^/2-1/2^ -^ enable us to conclude that 

(u^)^gN converges, up to an extraction, to some function u G L°°([0,T]; i?'^/^"^) fl 

Thus, we can pass to the limit in the system (3.7) and setting u := u + Uig. Then, 

(p, (7,u) satisfies the system (3.1). 

Step 5: Continuities in time. 

From the first equation of (3.1), we get dtp E L'^([0,T]; B^^"^"^'^^"^) which implies 
p G C([0,T]; i?'^/2~^''^/2). So does g in the same space. For u, we can derive, from 
the third equation of (3.1), that dtu G (L^ + L'^){[0,T]; B'i/2-i,rf/2) ^hj^h yields u G 
C([0,T];5'^/2-i'"'/2). 

3.4. Uniqueness. Let (pi, (7i, Ui) and (p2, g2, U2) be two solutions in F^ of (3.1) with 
the same initial data. Without loss of generality, we can assume that (p2,5'27U2) is 
the solution constructed in the previous subsection such that 

1 + inf P2it,x) > 0. 

(t,x)G[0,T]xR<^ 

We want to prove that (p2,fi'2, 112) = {pi,gi, ui) on [0,T] x M'^. To this goal, we shall 
estimate the discrepancy {6p,6g,6u) := (p2 — Pi,fi'2 — 91,^2 — ui) with respect to a 
suitable norm, satisfying 

' dtSp + U2 ■ V5p + 6u ■ Vpi = (5pdivu2 + (pi + l)div5u, 
dtSg + U2 ■ V6g + 6u ■ Vgi = —6gdivu2 — {gi + n)div5u, 
dt6u + U2 ■ V6u + 6u ■ Vui - (1 + p2)(pA(5u(p + A)Vdivu) (3.17) 

- (5p(pAui + (p + A)Vdivui) + Q{p2, 5-2) - Q{pi,gi) = 0, 
[{6p,6g,6u)\t=o = (0,0,0). 
We shall prove the uniqueness in a larger function space 

Ft:=(C([0,T];5^/2))1+i ^ (C{[0,T]; B''/^) f] L\[0,T]; B''/^+^))''. 
By Proposition 3.4, we get for all T' G [0, T] 

Pp|lL°°([0,r'];B''/2) 
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rT' 

^e^""'l'^^([o.T'i;s''/2+i) / [\\5u . VpiWsm + II 5pdivu2 11^^/2 

+ \\{pi + l)diY5u)\\Ba/2^dT 

<e^ll"2lLl([0.T'];fld/2+l) / ||5u||sd/2||pi 11^^/2+1 + || (5p || 5^/2 || U2 || B<*/2+l 

Jo ^ 

+ (1 + IIPi|lBd/2)||5u||5d/2+i dr. 

Using the Gronwall inequality, it yields 

PPllL°°([0,r'];_Brf/2) 

<e^ll"2ll.i([o,T'i;s^/2+i) / [||<5u||5./2||pi||5rf/2+i + (1 + ||pi||5./2)||5u||^./2+i]dr (3.18) 
Jo 

^Ct (||5u||i2([o^T'];Brf/2) + |l^u||ii([o,T'];_Bd/2+i)) , 

where Ct is independent of T'. 
Similarly, we have 

ll^fi'llL^([0,T'];Bd/2) ^ Ct (||5u||^2([o,T'];_Bd/2) + \\^^\\L^{[0,T'];Bd/2+i)) ■ (3.19) 

Applying Lemmas 3.2 and A. 3 to the third equation of (3.17), it yields 

II'^^IIl°°{[0,T'];B'*/2-1) + ll^u||^i([o,T'];B'*/2+i) 

rp/ 

+ ||Q(P2,^2) - <5(Pl,^l)||Bd/2-ijrfr. 

By Lemma A.l, we get 

\\Q{P2,92) -Q{Pi^9i)\\b"/2-i < (1+ ||(Pl,P2,^l,^2)||Bd/2)^(||5p||Bd/2 + ||5^||Bd/2)- 

Thus, it follows that 



(3.20) 



II^u|Il^([0,T'];B'^/2-1) + Pu||^i([o,T'];Brf/2+i) 
i^CriT' + r'^'^^)(||5p||ioo([o,T'];B'«/2) + P5'IIl°°([0,T'];B'*/2)), 

since ui e ^^([0, T]; 5^/^+1) ^y Lemma A.5. 

From (3.18)-(3.20), it yields, with the help of Lemma A. 5, that 

II^u|Il°°{[0,T'];B'^/2-1) + ll^u||^i([o,T'];B'*/2+l) 

^Ct(T' + T' )(||5u||i2([o,2./].Bd/2) + ||5u||ii([o^y/].5d/2+i)) 
^Ct(T' + T' ^)(||5u||^^([Q^y,j.^d/2-i) + ||5u||^i([o^j.,].^d/2+i)). 

Therefore, if we choose T' so small that Ct(T' + T'^/^) < 1, then we obtain that 
{6p,6g,6u) = (0,0,0) on the time interval [0,T']. As in the proof of uniqueness for 
global well-posedness, we can extend T' to T by the translation with respect to the 
time variable, i.e. (5p, 6g, 6u) = (0, 0, 0) on the time interval [0, T]. 
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3.5. A continuation criterion. 

Proposition 3.6. Under the hypotheses of Theorem 3.1, assume that the system (3.1) 
has a solution (p, (7, u) on [0, T) x W^ which belongs to F^, for all T' <T and satisfies 

P,^GL°°([0,T);5'^/2''^/2+1) jj^f p(t,x)>-l, f llVulloorft < 00. 

(t,x)G[0,T)xR'* Jo 

Then, there exists some T* > T such that {p,g,u.) may be continued on [0,T*] x R"' 
to a solution of (3.1) which belongs to F^,. 

Proof. Recall that u satisfies 

Ui + u ■ Vu - (1 + p)(pAu + (/i + A) Vdivu) + Q{p, g) = 0, u|(=o = "o- 

By Lemma 3.2, we get, for T' < T, that 

^r^ C-Jo^'f ||Vu||oc + ||p|P., ,,i)<it / II II , /" II II J4 

^6e ■'O V '^"fld/2+iy ||Uo||Bd/2-l,d/2 + / \\p\\Bd/2,d/2+ldt 

for some constant C depending only on d and viscosity coefficients. Thus, there exists 
a constant e > such that (3.1) with initial data {p{T — e),g{T — e), u(T — e)) yields 
a solution on [0,2£:]. Since the solution (p, 5f,u) is unique on [0,T), this provides a 
continuation of (p, g, u) beyond T. D 

Appendix A. Littlewood-Paley theory and Besov spaces 

This section is devoted to recall some properties of Littlewood-Paley theory and 
Besov spaces which will be used in this paper. For more details, one can see [6, 12] 
and references therein. 

Let ip : M.'^ ^ [0, 1] be a radial smooth cut-off function valued in [0, 1] such that 

r 1, lei ^ 3/4, 

^(0 = < smooth, 3/4 < |^| < 4/3, 
i 0, lel > 4/3. 

Let ip{^) be the function 

y.(e):=^(e/2)-^(e). 

Thus, ip is supported in the ball {.^ G M'' : |.^| ^ 4/3}, and ip is also a smooth cut-off 
function valued in [0, 1] and supported in the annulus {.^ G M"^ : 3/4 ^ |,^| ^ 8/3}. By 
construction, we have 

5^^.(2-^0 = 1, ve^o. 



One can define the dyadic blocks as follows. For /c G Z, let 

A,/ := ^-V(2~'0^/- 
The formal decomposition 

/ = $^A,/ (Al) 

fegZ 

is called homogeneous Littlewood-Paley decomposition. Nevertheless, (A.l) is true 
modulo polynomials, in other words (cf.[16]), if / G ^'(M.'^), then Ylkez ^kf converges 
modulo ,^[M'^] and (A.l) holds in y'{W^)/0^^% 
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Definition A.l. Let s G M. For / G ^'(^"')> we write 

fcez 

A difficulty comes from the choice of homogeneous spaces at this point. Indeed, 
II -Wb" cannot be a norm on {/ G ^'{M.'^) : ||/||bs < c>o} because ||/||_bs = means that 
/ is a polynomial. This enforces us to adopt the following definition for homogeneous 
Besov spaces (cf. [6]). 

Definition A. 2. Let s G M and m = -[d/2 + 1 - s]. If m < 0, then we define B'{R'^) 
as 

B-^ = {/ G y'{R'') : UWbs < oo and / = J^ A^/ in ^'(M^)}. 

fcez 

If m ^ 0, we denote by ^^ the set of d variables polynomials of degree less than or 
equal to m and define 

B' = [fE y\R^)l^ra ■■ WnW < OO and / = J] A,/ in J^'(M^)/^^}. 

fcgZ 

For the composition of functions, we have the following estimates. 

Lemma A.l. Let s > and u E B"^ D L°°. Then, it holds 

(i) Let F G PF/„'J+^'°°(M'^) with F(0) = 0. Then F{u) G B' . Moreover, there exists 
a function of one variable Cq depending only on s and F , and such that 

\\F{u)\\bs ^Cq{\\u\\loo)\\u\\b^. 

(ii) Ifu, V G B"^/^, iy-u) G B' for s E {-d/2, d/2] andG E ly/f ^1+^'°°(M'^) satisfies 
G"(0) = 0, then G{v) — G{u) E B^ and there exists a function of two variables C 
depending only on s, N and G, and such that 

\\G{v) - G{u)\\bs ^ C'dlwllioo, ||t;||L°=) (||M||Bd/2 + IK'llijd/2) \\v - u\\bs. 

We also need hybrid Besov spaces for which regularity assumptions are different in 
low frequencies and high frequencies [6] . We are going to recall the definition of these 
new spaces and some of their main properties. 

Definition A. 3. Let s, t eR. We define 

||/||bm=5^2^^||A,/||2 + 5^2^*||A,/||2. 

fc^O fc>0 

Let m = —[d/2 + 1 — s], we then define 

B''\R'^) ={f E ^'(M'^) : WfWBs.t < 00}, if m < 0, 

B''\R'') ={f E y'(R'')/^m : WfWBs.^ < 00}, if m ^ 0. 

Lemma A. 2. We have the following inclusions for hybrid Besov spaces. 
(i) We have B''' = B' . 

(ii) Ifs^t then B''^ = B' n BK Otherwise, B''^ = B' + BK 
(iii) The space B^'^ coincides with the usual inhomogeneous Besov space 5| ^ . 
(iv) Ifsi ^ S2 and ti ^ ta, then B'^^^^ -4 B'^'^\ 

Let us now recall some useful estimates for the product in hybrid Besov spaces. 
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Lemma A.3. Let si, S2 > and f,geL°°r] B'^''^ Then fg e B"^'"^ and 

Let s e (-C//2, d/2], f G B'^/'^ and g G B' , then fg G B' and 

llfgh^ < ||/||i?d/2||c/||B.. 

Let Si, S2, ti, t2 ^ d/2 such that min(si + S2, ti + t2) > 0, / G 5^^'*^ and g G -B^2-i2 ^ 
Then fg G ^^i+^^-i.ti+fe-i ^^^ 

||/5'llB»l+=2-d/2,ti+i2-d/2 < ||/||B=l.«l||fi'||B'=2.«2- 

For a, /3 G M, let us define the following characteristic function on Z: 

^«,/3., ^ _ f «, if fc ^ 0, 

Then, we can recall the following lemma. 

Lemma A. 4. Let F be an homogeneous smooth function of degree m. Suppose that 
—d/2 < Si,ti, S2,t2 ^ 1 + d/2. The following two estimates hold: 

\iFiD)Aki^r-Va),F{D)Aka)\ 

<7,2-'=(^"'"W-"^)||v||^./.+.||a||B^„.2||F(D)Afca||2, 
|(F(Z})Afc(v ■ Va), A,,6) + (Afe(v ■ V6), F(Z})A,a)| 

<7fc||v||B./2+i X {2-'^''''^''^\\F{D)A,ah\\b\\BH^^, 

+ 2"^-(^^^'^^W-™)||a||B=.^2||A,6||2), 

where (■, ■) denotes the 2-inner product, 'Ylike.zlk ^ 1 o-nd the operator F[D) is defined 
by F{D)f := ^-'F{0^f. 

In the context of this paper, we also need to use the interpolation spaces of hybrid 
Besov spaces together with a time space such as L^([0,T); 5*'*). Thus, we have to 
introduce the Chemin-Lerner type space (cf. [4]) which is a refinement of the space 
LP([0,T);5*'*). 

Definition A. 4. Let p G [1, oo], T G (0, oo] and Si, S2 G M. Then we define 

II/IIlp([0,T);B='*) = 2^2^1l^fc/||LP([0,T);L2) + / , 2^* || Afc/||^p([o^r).^2). 
fc^O fc>0 

Noting that Minkowski's inequality yields ||/||lp([o,t);_b»^*) ^ II/IIlp([o,t);B'=>«)' ^^ ^^~ 
fine spaces Lp([0,T); 5"'*) as follows 

ZP([0,r);5«'*) = {/ G L^([0,r);i?^'*) : ||/||i.([o,T);BM) < oo}. 

If T = oo, then we omit the subscript T from the notation L^([0,T); i?^'*), that is, 
L^(i?*'*) for simplicity. We will denote by C([0,T); S*'*) the subset of functions of 
L°°([0, T); S^'*) which are continuous on [0,T) with values in S'^'*. 

Let us observe that L\[0,T); B''^) = L\[0,T); B''^), but the embedding 

LP{[0,T);B''') cLP([0,r);5"'*) 

is strict if p > 1. 

We will use the following interpolation property which can be verified easily (cf. 

[1, 2]). 
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Lemma A. 5. Let s,t,si,ti,S2,t2 E^ and p,pi,p2 E [l,oo]. We have 



< 



1- 



where - = 77 + ^, s = 9si + {l- 9)s2 and t = 9ti + (1 - 9)12. 
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